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1. Introduction 

There has been considerable interest in studying quarkonia in hot medium since the 
pubhcation of the famous Matsui and Satz paper [1]. It has been argued that color 
screening in a deconfined QCD medium will suppress the existence of quarkonium 
states, signaling the formation of a quark-gluon plasma (QGP) in heavy-ion collisions. 
Although this idea was proposed a long time ago, first principle QCD calculations, which 
go beyond qualitative arguments, have been performed only recently. Such calculations 
include lattice QCD determinations of quarkonium correlators [21 |3l IH |5l E] ; potential 
model calculations of the quarkonium spectral functions with potentials based on lattice 
QCD [3 m m [101 [m [121 [13] (see also Ref. [H] for a review)), as well as effective 
field theory approaches that justify potential models and reveal new medium effects 
[T5l [T6l \TT\ UHl [19]. Spectral properties of heavy quark bound states are important 
ingredients in modeling of heavy quarkonium production in hot medium as will be 
discussed later. 

2. Color screening and deconfinement 

At high temperatures, strongly-interacting matter undergoes a deconfining transition 
to a quark-gluon plasma (QGP). This transition is triggered by a rapid increase of the 
energy and entropy densities, as well as the disappearance of hadronic states. (For 
recent reviews, see Ref. [20l [211 [22] ) . According to current lattice calculations, at zero 
net baryon density deconfinement occurs at Tc ~ (170 — 195) MeV [231 [21] Ijl 

J In the view of crossover nature of the deconfinement transition it is difficult to define the 
corresponding temperature interval precisely. Moreover due to discretization errors there is a 
discrepancy in the value of the deconfinement transition temperature obtained in calculations using 
different discretization schemes. Calculations with so-called stout staggered fermion action give a 
deconfinement transition temperature around 170MeV ^251 [26]. The new calculation with so-called 
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The QGP is characterized by color screening: the range of interaction between 

heavy quarks becomes inversely proportional to the temperature. Thus at sufficiently 

high temperatures, forming a bound state with a heavy quark (c or b) and its anti-quark 

becomes impossible. Color screening is studied on the lattice by calculating the spatial 
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Figure 1. Static quark singlet free energy versus quark separation calculated in 2+1 
flavor QCD on 16"^ x 4 lattices at different temperatures [36 (right). The combination 
r{Fi{r,T) ~ i^oo(T)) as function of rT (left). The sohd black line on the left plot is 
the parametrization of the zero temperature potential calculated in Ref. |24j . 
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Figure 2. Static quark free energy versus quark separation calculated in 2+1 flavor 
QCD on IG'^ X 4 lattices at different temperatures [35] (right). The combination 
r{F{r,T) - Foo{T)) as function of rT (left). The sohd black line on the left plot 
is the parametrization of the zero temperature potential calculated in Ref. |24) . 

correlation function of a static quark and anti-quark, which propagates in Euclidean time 
from T = to T = 1/T where T is the temperature (see Ref. [28] for a recent review). 
Two types of correlation functions are usually calculated on the lattice. The correlation 
function of Polyakov loops, which is also called the color averaged correlator 

G(r,r) = i(TrL(r)TrLt(0)), (1) 

where the temporal Wilson line is defined in terms of link variable Uo{xQ,r) as 
L(r) = n;^'=o^ Uo{xo,r), and the color singlet correlator 

Gi(r,T) = ^(TrL(r)Lt(0)) (2) 

HISQ action [57] seem to support this scenario. 
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The later is defined in Coulomb gauge [221 EOl ED E2] or in a gauge invariant but path 
dependent manner [28]. The singlet correlators recently have been also calculated in 
perturbation theory [3l]. The free energy of static quark anti-quark pair is defined 
through the logarithm of the color averaged correlator F{r, T) = — T In G{r, T). 
Analogously, the so-called singlet free energy is defined as Fi{r,T)/T = — In Gi (r, T) 
The singlet free energy in 2-1-1 flavor QCD is shown in Fig. [1] (right). The numerical 
results were obtained on 16'^ x 4 lattice using p4 action and quark masses which 
correspond to pion mass of about 220MeV [36]. As expected, in the zero temperature 
limit the singlet free energy coincides with the zero temperature potential calculated on 
the lattice [21]. Figure [H also illustrates that, at sufficiently short distances, the singlet 
free energy is temperature independent and equal to the zero temperature potential. At 
large distances it approaches a constant value Foo(T), which is twice the free energy of an 
isolated static quark. The range of interaction decreases with increasing temperatures. 
For temperatures above the transition temperature, T^, the heavy quark interaction 
range becomes comparable to the charmonium radius. In particular, the combination 
r(Fi(r, T) — Foo(T)) is expected to decay exponentially at large distances, rrriD > 1, 
with rriD being the Debye mass. Indeed, for distances rT > 0.8fm the screening is 
exponential as can be seen from Fig. [1] (right). Based on this general observation, one 
would expect that the charmonium states, as well as the excited bottomonium states, 
do not exist above the deconfinement transition. (In the literature, this is often referred 
to as dissociation or melting). 

The path/gauge dependence of singlet free energy makes its interpretation 
somewhat complicated. Therefore, on the lattice one also calculates the Polyakov loop 
correlator. The results of the calculations in 2-1-1 flavor QCD with p4 action on 16'^ x 4 
lattices are shown in Fig. |2]in terms of the logarithm of the correlation function, i.e. the 
physical free energy of static quark anti-quark pair. Compared to the singlet free energy 
the physical free energy shows much stronger temperature dependence and is different 
from the zero temperature potential at all temperatures. In the low temperature region 
the strong temperature dependence of the free energy can be understood as resulting 
from significant contribution from higher excited states [33]. At high temperatures, 
the origin strong temperature dependence is more difficult to understand. In the short 
distance limit within the effective field theory framework it can be attributed to the 
contribution from octet degrees of freedom [35] . As one can also see from Fig. |2] the free 
energy of static QQ pair also shows exponential screening at large distances. Thus the 
strong color screening is a real effect. In summary, lattice calculations of static quark 
anti-quark correlators provide evidence that in the deconfined phase there are strong 
screening effect at distances relevant for quarkonium physics. In the next section the 
implication of these findings on quarkonium spectral functions will be discussed. 

§ The identification of the correlators defined by Eqs. (H]) and ^ with color averaged and color singlet 
channels is based on perturbative arguments, e.g. see discussion in Ref. |37| . Non-perturbativcly such 
identification becomes problematic as was shown in Ref. p8]. At short distances effective field theory 
concepts can be invoked to define color singlet and color averaged channels [35) . 



Quarkonium in Hot Medium 



4 



3. Quarkonium spectral functions 

In-medium quarkonium properties are encoded in the corresponding spectral functions, 
as are their dissolution at high temperatures. Spectral functions are defined as the 
imaginary part of the retarded correlation function of quarkonium operators. Bound 
states appear as peaks in the spectral functions. The peaks broaden and eventually 
disappear with increasing temperature. The disappearance of a peak signals the melting 
of the given quarkonium state. 

In lattice QCD, the meson correlation functions, G{t,T), are calculated in 
Euclidean time. These correlation functions are related to the spectral functions cr{u, T) 
as 



Detailed information on G{t, T) could make it possible to reconstruct the spectral 
function from the lattice data. In practice, however, this turns out to be very difficult 
task because the time extent is limited to 1/T, see the discussion in Ref. [5] and references 
therein. 

The quarkonium spectral functions can be calculated in potential models using the 
singlet free energy from Fig. [1] or with different lattice-based potentials obtained using 
the singlet free energy as an input [I2ill3j. The results of calculations in quenched QCD 
are shown in Fig. [3] for S'-wave charmonium and bottomonium spectral functions |12j . 
All charmonium states are dissolved in the deconfined phase, while the bottomonium 
IS state may persist up to temperature of about 2Tc. The temperature dependence 
of the Euclidean correlators can be predicted using Eq. ([3]) and the calculated spectral 
functions. On the lattice the temperature dependence of the correlation function is 
studied in terms of the ratio G{T,T)/Grec{T,T), where 



If the spectral function is temperature independent G/Grec = 1- It turns out, somewhat 
surprisingly, that the Euclidean correlation functions in the pseudo-scalar channel show 
very little temperature dependence irrespective of whether a state remains bound (the 
r]b{lS) )or not (the ?7c(l>S')). Note also that correlators from potential models are in 
accord with the lattice calculations (see insets in Fig. |3]). Initially, the weak temperature 
dependence of the correlators was considered to be evidence for the survival of different 
quarkonium states [4]. It is now clear that this conclusion was premature. There is a 
large enhancement in the threshold region of the spectral functions relative to the free 
spectral function, as shown in Fig. |3l This threshold enhancement compensates for the 
absence of bound states and leads to Euclidean correlation functions with a very weak 
temperature dependence [12]. It further indicates strong residual correlations between 
the quark and anti-quark, even in absence of bound states. Present lattice calculations 
of the spectral functions cannot discriminate between bound state peaks and threshold 
enhancement. 
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Figure 3. The 5-wave charmonium (left) and bottomonium (right) spectral functions 
calculated in potential models |12j. Insets: correlators compared to lattice data [41) . 
The dotted curves are the free spectral functions. 



Similar analysis has been done for the P-wave charmonium and bottomonium 
spectral functions [121 US]- Here the contribution of the zero mode is important. The 
zero mode contribution is responsible for large enhancement and strong temperature 
dependence of G/Grec in the scalar and axial- vector channels [121 ES]- A smaller 
temperature dependence of G/Grec due to zero mode was observed in the vector channels 
as well |l0]. Here the zero mode contribution is related to heavy quark diffusion. A 
comprehensive study of the zero mode contribution was presented in Ref. [4Ij where 
it was concluded that almost the entire temperature dependence of the correlation 
function is due to this contribution. Furthermore, the zero mode contributions which 
are related to generalized susceptibilities are well described by quasi-particle model 
[41] . This provides additional input for in- medium properties of heavy quarks, e.g. the 
temperature dependent heavy quark mass, which will be important for constructing 
more realistic potential model at finite temperature. While the above studies have been 
performed in quenched approximation, recently attempts to extend the quasi-particle 
model in the charm quark sector to full QCD have been presented in Ref. [12]. 

To estimate the dissociation temperatures, i.e. the temperature above which no 
resonance like peak are observed in the quarkonium spectral functions, it is not sufficient 
to consider only the effects due to color screening. One has to consider also the effect of 
thermal broadening of different states. Upper bound on the dissociation temperatures 
can be obtained from the analysis of the spectral functions calculated in potential model 
and using simple phenomenological estimates for the thermal width [13j. This analysis 
shows that most of quarkonium states dissolve at temperatures < 1.2Tc and only the 
ground state bottomonium may survive up to temperatures as high as 2Tc. 

The application of potential models can be justified using an effective field theory 
approach. The existence of distinct energy scales related to the heavy quark mass m, 
the inverse size mv (where v is the heavy quark velocity), and the binding energy mv'^ 
makes it possible to construct a sequence of effective theories at zero temperature. The 



Quarkonium in Hot Medium 



6 



effective theory which emerges after integrating out the scales m and mv is pNRQCD. 
The Lagrangian of this effective theory contains singlet and octet meson fields composed 
of heavy quarks, which are coupled to soft gluon fields at scale mv"^. pNRQCD is 
equivalent to the potential model at T = at leading order, where the coupling to 
the octet channel can be neglected. It is possible to extend this approach to finite 
temperature where additional scales, the temperature T, the Debye mass m^) ~ gT, 
and the magnetic scale g'^T are present. In the weak coupling regime, (7^1, these 
scales are well separated. Depending on how the thermal scales are related to the zero 
temperature scales, the various hierarchies make it possible to derive different effective 
theories for quarkonium bound states at finite temperature [I9]. In the weak-coupling 
QCD approach, thermal corrections to the potential are obtained when the temperature 
is larger than the binding energy. An important general result of these effective theories 
is that the potential acquires an imaginary part [IHl HEl HTJ [TBI HQ] . The imaginary part 
of the potential smears out the bound state peaks of the quarkonium spectral function, 
leading to their dissolution prior to the onset of Debye-screening in the real part of the 
potential (see e.g. the discussion in Ref. [18]). 

4. Dynamical models for quarkonium production 

While knowing the quarkonium spectral functions in equilibrium QCD is necessary, 
it is insufficient to predict effects on their production in heavy-ion collisions because, 
unlike the light degrees of freedom, heavy quarks are not fully thermalized in heavy-ion 
collisions. Therefore, it is nontrivial to relate the finite temperature quarkonium spectral 
functions to quarkonium production rates in heavy ion collisions without further model 
assumptions. The bridge between the two is provided by dynamical models of the 
matter produced in heavy-ion collisions. Some of the simple models currently available 
are based on statistical recombination |13]; statistical recombination and dissociation 
[HIH5] : or sequential melting |16]. Here we highlight a more recent model, which makes 
closer contact with both QCD and experimental observations [17]. 

The bulk evolution of the matter produced in heavy-ion collisions is well modeled 
by hydrodynamics, see Ref. [IBJ for a recent review. The large heavy quark mass makes 
it possible to model its interaction with the medium by Langevin dynamics [19]. Such 
an approach successfully describes the anisotropic flow of charm quarks observed at 
RHIC [m [50l [51] (see also Ref. [52] for a recent review). Potential models have shown 
that, in the absence of bound states, the QQ pairs are correlated in space [121 |T3]. This 
correlation can be modeled classically using Langevin dynamics, including a drag force 
and a random force between the Q (or Q) and the medium as well as the forces between 
the Q and Q described by the potential. It has been shown that a model combining 
an ideal hydrodynamic expansion of the medium with a description of the correlated 
QQ pair dynamics by the Langevin equation may explain why, despite the fact that 
a deconfined medium is created at RHIC, there is only a 40 — 50% suppression in the 
charmonium yield. The attractive potential and the finite lifetime of the system prevents 
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the complete decorrelation of some of the QQ pairs [H]. Once the system has sufficiently 
cooled, these residual correlations make it possible for the Q and Q to form a bound 
state. The details of the suppression pattern depend on the model parameters, such as 
Tc and the way recombination is implemented. The later appears to be important in 
the central collisions |53]. However, the results are not very sensitive to the choice of 
the potential [M] . 

The above approach, which neglects quantum effects, is applicable only if there are 
no bound states, as it is likely to be the case for the J/ip. If heavy quark bound states are 
present, as is probable for the T{1S), the thermal dissociation rate will be most relevant 
for understanding the quarkonium yield. It is expected that the interaction of a color 
singlet quarkonium state with the medium is much smaller than that of heavy quarks. 
Thus, to first approximation, medium effects will only lead to quarkonium dissociation. 

5. Concluding remarks 

Lattice calculations of static quark anti-quark correlators provide evidence for strong 
screening effects in the deconfined phase. Potential model calculations based on lattice 
QCD, as well as re-summed perturbative QCD calculations indicate that all charmonium 
states and excited bottomonium states dissolve in the deconfined medium. This leads 
to the reduction of the quarkonium production yield in heavy-ion collisions compared to 
binary-scaling of pp collisions. Survival of residual correlations in the deconfined phase 
as well as recombination effects at lower temperatures lead to non-zero quarkonium yield 
in heavy ion collisions. 

It turns out, that despite melting of quarkonium states corresponding meson 
correlation function in Euclidean time show very little temperature dependence. Almost 
the entire temperature dependence of the Euclidean correlation functions is due to the 
zero mode contribution. Statements about the existence of heavy quark bound states 
based on lattice correlation functions and spectral functions need to be revisited. Yet, 
lattice results on correlation functions are still valuable tools for constraining potential 
models. 

One of the great opportunities of the LHC and RHIC-II heavy-ion programs is 
measurement of bottomonium yields. From a theoretical perspective, bottomonium 
is an important and clean probe for at least two reasons. First, the effective field 
theory approach, which provides a link to first principle QCD, is more applicable for 
bottomonium due to better separation of scales and higher dissociation temperatures. 
Second, the heavier bottom quark mass reduces the importance of recombination effects, 
making bottomonium a good probe of dynamical models. 
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